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THE  USE  OF  SATELLITE-OBSERVED  CLOUD  PATTERNS  IN 
NORTHERN  HEMISPHERE  500-MB  NUMERICAL  ANALYSIS 


Roland  E.  Nagle  and  Christopher  M.  Hayden 


ABSTRACT.   A  quasi-objective  method  for  deriving  500-irib 
geopotential  heights  with  the  aid  of  satellite  cloud 
observations  is  presented.  The  method  uses  satellite- 
observed  cloud  patterns  in  conjunction  with  a  forecast 
500-mb  height  field  which  is  separated  into  additive 
short-  and  long-wave  component  fields.   Empirical 
relationships  between  the  cloud  patterns  and  the  short- 
wave component  field  are  used  to  modify  the  pattern  of 
the  latter,  and  regression  equations  are  used  to  specify 
extrema  values.  The  500-mb  heights  are  retrieved  by  the 
direct  addition  of  the  modified  short-wavelength  field  to 
the  long -wavelength  field.  Procedures  for  using  this 
method  operationally  are  presented,  and  results  of  real- 
time application  are  shown.  It  is  concluded  that  the 
resulting  modifications  can  contribute  to  significant 
improvements  in  the  analysis. 

1  .   INTRODUCTION 

Prerequisite  to  the  use  of  satellite  cloud  data  in  numerical  weather 
analysis  is  the  establishment  of  consistent  relationships  between  cloud 
patterns  and  parameters  which  can  be  used  in  or  readily  transformed  into 
data  suitable  for  a  numerical  analysis  system.   The  literature  is  abundant 
with  reports  on  the  relationships  between  cloud  pattern  features  and 
descriptive  features  of  synoptic  meteorology,  i.e.,  fronts,  troughs, 
ridges,  jet  streams,  etc.  The  results  of  these  studies,  which  were  first 
consolidated  by  Widger  (I96l4.)j  provide  a  useful  background  for  inter- 
preting satellite-observed  cloud  patterns,  but  do  not  offer  a  means  for 
utilizing  these  data  in  numerical  weather  analysis. 

An  approach  to  this  problem  was  first  proposed  by  Bristor  and 
Ruzecki  (1960).   They  suggest  that  use  be  made  of  the  approximate  rela- 
tionship between  clouds  and  vertical  motion  at  500-mb  as  expressed  in 
the  vorticity  advection  term  of  the  quasi-geostrophic  "Omega  Equation." 
The  procedure  required  the  modification  of  the  height  or  the  stream 
function  field  (or  their  Laplacians)  at  the  500-mb  level  on  the  basis  of 
satellite-observed  cloud  patterns.  McClain,  Brodrick,  and  Ruzecki  (1965? 
1966)  attempted  reanalyses  using  these  relationships,  and  Bradley,  Hayden, 
and  Wiin-Nielsen  (I966)  used  this  approach  in  further  experiments. 


The  approach  as  reported  by  these  investigators  has  certain  limita- 
tions .  Although  modifying  the  Laplacian  field  changes  the  vorticity 
advection,  the  effects  are  not  known  until  after  the  height  field  or 
stream  function  field  is  retrieved  by  relaxation.   Relaxation  also  ob- 
structs preservation  of  conventional  data  and,  if  total  vorticity  is  not 
conserved  within  the  reanalysis  area,  can  cause  changes  outside  the  re- 
analysis  area  when  the  height  or  stream  function  field  is  retrieved.   The 
magnitude  of  modifications  in  the  Laplacian  field  is  critical  and  cannot 
be  established  objectively.  Finally,  care  must  be  taken  in  the  reanalysis 
to  insure  that  negative  absolute  vorticities  are  not  produced. 

The  approach  presented  in  this  report  is  an  extension  of  Bristor's 
and  Ruzecki's  (I96O)  suggestions  with  several  significant  variations. 
Also,  it  is  a  simplification  of  a  system  first  proposed  by  Nagle,  et.  al., 
(1966)  and  developed  by  Nagle  and  Clark  (1968).   The  approach  is  analogous 
to  perturbation  theory  in  that  it  requires  the  objective  separation  of  the 
500-mb  height  field  into  additive  long-  and  short-wavelength  components, 
i.  e.,  the  spatial  mean  flow  and  superimposed  disturbances.   The  good 
relationship  found  between  the  cloud  patterns  and  certain  features  of  the 
short-wavelength  field  permits  manual  modification  of  the  latter  in  areas 
of  sparse  conventional  data.   The  use  of  relaxation  techniques  to  retrieve 
the  500-mb  height  values  is  not  required  because  the  direct  addition  of 
the  short-  and  long -wavelength  components  produces  the  modified  500-mb 
height  field.   Statistical  methods  are  used  to  derive  regression  equations 
from  which  the  magnitude  of  extrema  in  the  short-wavelength  field  can  be 
computed  objectively.   The  modifications  effected  in  the  total  height 
field  are  relatively  minor  and  do  not  significantly  change  the  pattern  of 
the  long -wavelength  field.   Thus,  the  vorticity  advection  at  500  mb  is 
also  adjusted  to  conform  to  the  observed  cloud  patterns  (Nagle  and  Clark, 
1968).   The  results  indicate  that,  with  proper  interpretation  of  the  cloud 
patterns,  the  short -wavelength  features  of  the  500-mb  height  field  can  be 
correctly  introduced  into  numerical  analysis. 

2.   SEPARATION  BY  SCALE 

The  method  for  objectively  separating  meteorological  fields  into 
additive  components  was  developed  by  Holl  (I963).   The  technique^ essen- 
tially an  objective  version  of  Fjortoft's  (19^2)  method,  successively 
applies  a  smoothing  operator  to  a  numerical  field  until  the  amplitude  of 
a  specified  wavelength  component  is  reduced  to  a  specified  percentage  of 
its  initial  value.   Since  the  operator  is  linear,  all  additions  and 
subtractions  to  yield  components  of  the  initial  field  are  commutative. 

The  smoothing  is  defined  by  (Holl,  I963): 

ZlW  =  Zq  +  cT^V^Z  da  =  Zq  -  Zsw 
»/0 

CD 


where  Z-qt  is  the  long-wavelength  field,  Zq  is  the  original  field,  Zg^  is 
the  short-wavelength  field,  C  is  a  constant,  and  a  is  a  parameter  repre- 
senting the  degree  of  smoothing.  Appendix  B  contains  a  detailed  discus- 
sion of  the  smoothing  procedure;  Appendix  C  is  the  computer  program  which 
perfonns  the  smoothing. 

The  degree  of  smoothing  to  be  applied  is  somewhat  arbitrary,  but  there 
is  evidence^  that  a  reasonable  separation  can  be  made  between  the  short- 
and  long-wave  components  at  500  mb  with  the  degree  of  smoothing  set  at 

ct  =  5.  For  this  amount  of  smoothing  the  amplitude  of  a  system  with  a 
wavelength  of  7  gridlengths  (1  gridlength  =  38I  km  at  60°  N.)  is  reduced 
to  5  percent  of  its  initial  value.  A  system  of  this  size  corresponds 
approximately  to  zonal  wave  number  10  at  h^°   N.   In  contrast,  the  ampli- 
tude of  a  system  with  wave  number  3  (23  gridlengths)  is  reduced  by  less 
than  k   percent.  The  degree  of  smoothing  used  throughout  this  study  was 

a  =  ^.      The  corresponding  amplitude  reduction  as  a  function  of  wavelength 
(in  gridlengths)  is  shown  in  figure  1 . 

Figures  2  through  h   are  examples  of  the  fields  produced  by  the  scale 
separation  technique.   Figure  2  is  an  initial  500  mb  height  field;  figure 
3  is  the  field  produced  by  the  smoothing  process,  i.  e.,  the  long -wave 
pattern.  By  subtracting  the  long -wave  field  from  the  initial  500-mb 
analysis  at  gridpoints,  the  components  which  have  been  smoothed  out  of 
the  initial  analysis  can  be  retrieved.  The  short-wavelength  component 
field  corresponding  to  figures  2  and  3  is  shown  in  figure  U. 

3.   SIGNIFICANCE  OF  THE  SHORT -WAVELENGTH  COMPONENT 

The  short-wavelength  component  field  is  significant  because  it  approx- 
imates the  field  of  relative  vorticity.   Certain  characteristic  features 
of  the  cloud  patterns  are  related  to  the  500-mb  relative  vorticity  field 
and  therefore  to  the  short-wavelength  field  (see,  for  example,  Nagle, 
et.  al.,  1966).   Thus,  configurations  of  the  satellite-observed  cloud 
field  may  be  used  to  modify  the  short-wavelength  component  field.  These 
short-wave  modifications  are  in  turn  easily  incorporated  in  the  total 
height  field  by  simple  addition  to  the  long-wavelength  component  field. 
We  shall  first  discuss  the  relationship  between  the  short-wave  and 
relative  vorticity  fields  and  then  describe  relationships  between  the 
relative  vorticity  and  the  cloud  patterns. 

From  (1)  the  short-wavelength  component  is  given  by 

2sw  =-Cr°72zda 

^0  (2) 


1  Meteorology  International  Inc.,  Quarterly  Performance  Report  No.  h, 
Contract  No.  N00228-67-C-2759,  Monterey,  Calif.,  September,  1 968 . 
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Figure  1. — Filter  curve  for  smoothing  (a  =  ^) 
One  gridlength  equals  36 1  km  at  60°  N. 


Figure   2.— 500-mbgeopotential  height   field    (decameters),    1200  GMT, 
March  l6 ,   1970. 


Figure   3. ~500-mb long-wavelength   field    (decameters),   1200  GMT. 
March  l6 ,   1970. 
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Figure  k. —  500-mb short-wavelength  field  (meters),  1200  GMT, 
March  l6 ,  1970. 


-120 


H+ 


\ 


.+60 


N 


L- 

-120 


/ 


/ 


-60' 


y 


y 


>v. 


H+ 


~+60^ 


Figure  5» — 500-mb short-wavelength  field  (meters),  1200  GMT, 
February  8,  1968. 


figure  6. — 5 00 -mb relative  vorticity  field  (xlO"-'  sec"-'-),  1200  GMT, 
February  8,  1968. 
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By  the  mean  value  theorem  for  integrals; 


(3) 


where  Z(a  -/-)   represents  the  partially  smoothed  field  at  some  degree  of 
smoothing,  a   -",  between  0  and  a   .   The  geos trophic  relative  vorticity  q 
is  given  by 


q  =  gf-^  V^  Zq 


(U) 


where  g  is  the  acceleration  of  gravity  and  f  is  the  coriolis  parameter. 
Therefore,  the  correspondence  between  the  relative  vorticity  and  the 
short-wavelength  component  depends  on  the  latitude  and  the  equivalence  of 

Because  the  del  square  operator  is  wavelength  dependent,  there  is  no 
reason  to  expect  direct  proportionality  between  the  magnitudes  of  the 
fields  V  Z  (ct-x-)  and  V  Zq  •  Nevertheless,  the  patterns  of  these  fields 
should  be  very  similar.  In  particular,  the  locations  where  the  fields 
have  either  extrema  or  zero  values  should  be  similar.  Figures  $   and  6  are 
the  relative  vorticity  and  short-wavelength  component  fields  for  1200  GMT, 
February  8,  1968.  They  show  that  the  short-wave  field  retains,  with  light 
smoothing,  the  patterns  of  the  relative  vorticity  field.  Note  that  the 
short-wavelength  field  is  opposite  in  sign  from  the  vorticity  field. 

The  correspondence  found  between  identifiable  cloud  features  and 
synoptic  entities  associated  with  specific  details  of  the  $00-mb   relative 
vorticity  field  is  illustrated  in  figure  7  and  can  be  summarized  as 
follows : 

1 .  The  boundaries  of  cellular  convective  cloud  areas  and  the  trailing 
edges  of  frontal  cloud  bands  correspond  to  the  zero  lines  in  the  relative 
vorticity  field. 

2.  The  locations  of  spiral  cloud  centers  associated  with  cold-core 
systems  correspond  to  positive  vorticity  maxima. 


The  keying  in  figure  7  is  applicable  to  all  subsequent  figures  in  this 

report. 
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Figure  7. — Schematic  representation  of  the  relationships  between 
cloud  patterns  and  features  of  the  ^00-mb   relative 
vorticity  field. 
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3 .   The  leading  edges  of  frontal  cloud  bands  correspond  to  axes  of 
maximum  negative  vorticity. 

h'  The  points  marking  the  apex  of  the  anticyclonic  curvature  on  the 
leading  edges  of  frontal  cloud  bands  correspond  to  negative  vorticity 
maxima. 

These  relationships  permit  a  rather  complete  specification  of  the 
pattern  of  the  short -wavelength  field  from  an  interpretation  of  the  cloud 
patterns.  It  remains,  .however,  to  specify  the  magnitude  of  the  negative 
extrema  in  the  short-wavelength  field. 

h.      STATISTICAL  SPECIFICATION  OF  THE  CENTRAL  VALUES 
OF  SHORT-WAVELENGTH  SYSTEMS 

A  statistical  investigation  was  undertaken  to  determine  relationships 
between  the  magnitude  of  the  negative  maxima  in  the  short-wavelength  field 
and  parameters  derived  from  both  the  long-wavelength  field  and  measure- 
ments of  the  cloud  systems.   The  concept  was  to  derive  these  relationships 
from  a  dependent  sample  over  a  relatively  dense-data  area  so  that  the 
resulting  regression  equations  could  be  used  to  specify  the  magnitude  of 
systems  over  sparse  data  regions.   The  dependent  data  were  derived  over 
the  North  Atlantic  Ocean,  in  the  approximate  region  N.  25° -70°,  E.  10°- 
W.  80°. 

The  cloud  patterns  were  analyzed  as  follows  (fig-  8).   The  location 
of  a  spiral  cloud  center  was  entered  onto  a  1:30,000,000  polar  stereo- 
graphic  base  map.  A  curved  line  was  drawn  radially  outward  from  the 
center  along  the  leading  edge  of  the  outermost  cyclonically  curved  band 
(center  to  M' ) .   In  the  eastern  quadrant  of  the  cloud  pattern,  the  line 
was  continued  along  the  leading  edge  of  the  brightest-appearing  cloudi- 
ness and  merged  with  the  leading  edge  of  the  frontal  band.   This  line  was 
continued  along  the  frontal  band  to  the  point  where  the  500-mb  trough 
line  intersected  the  band  (M'  to  T),  i.e.,  to  the  point  where  the  frontal 
cloud  band  either  ends  or  changes  abruptly  in  character.   Next,  a  base 
line  was  drawn  from  the  center  to  the  point  T,  where  the  500-mb  trough 
intersected  the  cloud  band.   The  length  of  this  line  is  defined  as  the 
cloud  amplitude.  A  second  line  was  then  drawn  perpendicular  to  the  base 
line  at  a  location  in  which  the  distance  along  this  second  line  through 
the  leading  and  trailing  edges  of  the  curve  depicting  the  cloud  pattern 
was  maximized.   This  line  (M  to  M' )  is  defined  as  the  maximum  diameter 
of  the  cloud  system.  Measurements  of  both  cloud  amplitude  and  maximum 
diameter  were  made  in  tenths  of  gridlengths.   These  measurements  were 
obtained  for  each  cloud  system  that  could  be  associated  with  a  specific 
maximum  in  the  short-wave  field.   In  addition  to  the  cloud  parameters, 
the  following  were  computed  at  the  location  of  each  center:   the  latitude, 
the  height  value  and  the  Laplacian  of  the  long -wavelength  field,  the  u 
and  V  component  of  the  geostrophic  wind  (derived  from  the  long  wavelength 
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Figure  8. — Geometrical  construction  for  determining  the  diameter 
and  amplitude  of  a  cloud  system. 
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field),  and  the  long-wave  thickness  advection.  No  attempt  was  made  to 
classify  the  cloud  systems  by  stage  of  development. 

The  dependent  sample  consisted  of  1^6  cases  over  a  3-i7ionth  period 
CFebruary  and  October  1 968  and  June  I969).   The  population  of  cases  was 
analyzed  statistically  by  means  of  a  stepwise  linear  regression  screening 
procedure.  After  experimentation,  the  following  predictors  were  selected 
for  operational  usage:   latitude,  the  finite  difference  Laplacian  of  the 
long-wavelength  field,  and  the  maximum  diameter  and  amplitude  of  the  cloud 
system.  These  were  chosen  because  they  produced  the  greatest  reduction  of 
variance  among  the  parameters  which  are  most  amenable  for  manual,  opera- 
tional use.  A  second  regression  equation  with  the  cloud  parameters 
excluded  was  derived  from  this  sample.  This  equation  is  used  when  spec- 
ification of  the  cloud  parameters  is  ambiguous .   The  coefficients  of  the 
regression  equations  are  provided  in  table  1 . 

Table  1 . — Regression  equation  coefficients 
Variable        With  cloud  parameters   Without  cloud  parameters 

Constant +21.5  -^2.6 

Latitude -1.7  -0.8 

Laplacian -2.1  -2.5 

Max  diameter 10. I|. 

Amplitude -1.0,' 


5.   CONSTRUCTION  OF  500-MB  HEIGHTS 
AT  THE  LOCATION  OF  SPIRAL  CLOUD  CENTERS 

Since  parameters  derived  from  the  long-wavelength  field  are  used  in 
the  regression  equation  for  computing  the  intensity  of  the  maxima  in  the 
short-wavelength  field,  and  since  the  retrieval  of  absolute  500-mb  heights 
requires  the  heights  of  the  long-wavelength  field,  it  must  be  assumed  that 
the  long-wavelength  field  is  known.  In  a  realtime  operation  the  long- 
wavelength  field  derived  from  the  current  analysis  is  not  available  when 
the  satellite  observations  are  being  interpreted.  It  is  assumed  that 
there  is  sufficient  information  in  the  1 2 -hour  numerical  forecast  (the 
"first-guess"  field)  to  allow  the  accurate  specification  of  the  long- 
wavelength  pattern  in  the  500-mb  height  field.  This  assumption  has  been 
tested  by  determining  the  magnitude  of  the  differences  between  the  long- 
wavelength  fields  derived  from  first-guess  fields  and  operational 
analysis  fields . 
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The  heights  of  the  long -wavelength  field  as  derived  from  the  first- 
guess  field  and  as  derived  from  the  analyzed  fields  were  compared  daily 
over  a  3 -month  period.  The  same  comparison  was  made  for  the  Laplacian 
of  the  long -wavelength  field.   These  comparisons  were  made  over  the 
North  Atlantic.  For  each  day  in  this  sample,  the  following  parameters 
were  derived:   The  maximum  and  minimum  values,  the  mean  value,  and  the 
standard  deviation.  Also  the  correlation  coefficient,  the  absolute 
maximum  difference,  and  the  root-mean-square  (RMS)  error  between  the 
first-guess  field  values  and  the  analyzed  values  were  computed.  Tables 
2  and  3  are  statistics  for  each  of  the  three  months. 

By  applying  the  coefficients  of  table  1  to  the  corresponding  values 
of  tables  2  and  3  the  errors  involved  in  accepting  the  long-wavelength 
pattern  derived  from  the  first-guess  field  can  be  estimated.  Applying 
worst-on-worst  principles,  the  extreme  error  in  the  computed  500-rab 
heights  at  the  location  of  the  vorticity  centers  attributable  to  the 
errors  in  the  long -wavelength  field  should  be  no  greater  than  80  meters. 
But  considering  the  mean  RMS  errors,  a  more  reasonable  value  for  the 
expected  error  should  be  on  the  order  of  20-30  meters.  This  is  considered 
acceptably  low  for  use  in  the  current  context. 

By  accepting  the  long-wavelength  field  derived  from  the  first-guess 
field  as  a  reasonable  base,  the  500-mb  heights  are  retrieved  by  the 
direct  addition  of  the  two  scale  components. 

The  500-mb  heights  at  the  location  of  the  spiral  cloud  centers  are 
obtained  by  interpreting  the  cloud  patterns  in  terms  of  the  parameters 
required  in  the  solution  of  the  regression  equation  (fig.  8).   The 
Laplacian  of  the  long -wavelength  field  is  computed  at  the  grid  point 
closest  to  the  location  of  the  spiral  cloud  center.  The  central  value  of 
the  associated  short-wavelength  system  is  then  retrieved  by  solving  the 
regression  equation.  Adding  this  value  to  the  corresponding  height  in 
the  long  wavelength  field  produces  the  $00-mb   height  at  the  location  of 
the  spiral  cloud  center. 

Where  the  short-wavelength  field  has  a  zero  value  (the  zero  line), 
the  value  of  the  long -wavelength  is  the  500-mb  height.  An  example  of 
this  is  illustrated  in  figure  9.  As  shown  in  figure  9  (left),  the 
location  of  the  zero  line  in  the  short-wavelength  field  is  inferred 
from  the  cloud  patterns .   The  zero  line  is  superimposed  upon  the  long- 
wavelength  field,  and  500-mb  heights  can  be  read  at  each  intersection 
of  the  zero  line  with  the  contours  of  the  long -wavelength  field. 


16 


Table  2. — Statistics  from  comparisons  of  long -wavelength  heights 
as  derived  from  the  "first  guess"  field 
and  the  analyzed  $00-mb   height  field 


Absolute  value  (meters): 

Analysis 

Guess 

Standard  deviation  (meters) 

Analysis 

Guess 

Maximujn  difference  C^fieters)' 

Correlation  coefficient 

RMS  error  (meters) 


Absolute  value  (meters): 

Analysis 

Guess 

Standard  deviation  (meters) 

Analysis 

Guess 

Maximum  difference  (meters)- 

Correlation  coefficient 

RMS  error  (meters) 


Absolute  value  (meters): 

Analysis 

Guess 

Standard  deviation  (meters) 

Analysis 

Guess 

Maximum  difference  (meters)- 

Correlation  coefficient 

RMS  error  (meters) 


Extreme 

Monthly 

Extreme 

maximum 

mean 

minunim 

February  1 968 

592U 

5669 

i|836 

5929 

5677 

i;8U8 

307.3 

2I4I  .1 

175.0 

296.8 

236.0 

168.1; 

6U.8 

h  .0 

23.9 

.999 

.998 

.996 

22. U 

16.2 

6.6 

October  1968 

$933 

5720 

5161 

59i|1 

5723 

5155 

215.5 

178.3 

127.7 

211  .h 

177.5 

.    128.I; 

39.3 

25.0 

lU.l 

.996 

.989 

.975 

15.0 

10.3 

5.6 

February  1969 

5896 

5616 

l;9l;1 

5900 

5618 

1;951 

2li7.1 

201.0 

150.0 

239.0 

197.6 

11;8.1; 

5I1.7 

37.3 

21;. 6 

.999 

.998 

.996 

21.5 

11;. 8 

7.9 
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Table  3. — Statistics  from  comparisons  of  the  Laplacian  of  the  lon^ 
wavelength  height  field  as  derived  from  the  "first  guess"  field 
and  the  analyzed  500-mb  height  field 


Absolute  value: 

Analysis 

Guess 

Standard  deviation: 

Analysis 

Guess 

Maximum  difference 

Correlation  coefficient- 
RMS  error 


Absolute  value: 

Analysis 

Guess 

Standard  deviation: 

Analysis 

Guess 

Maximum  difference 

Correlation  coefficient- 
RMS  error 


Absolute  value: 

Analysis 

Guess 

Standard  deviation 

Analysis --- 

Guess 

Maximum  difference 

Correlation  coefficient- 
RMS  error 


Extreme 

Monthly 

Extreme 

maximum 

mean 

minunum 

February  I968 

+85.0 

+0.6 

-50.9 

+8U.0 

+1.9 

-50.1 

26.2 

20.9 

13.2 

26.1 

20.5 

12.9 

11.5 

8.2 

5.2 

.998 

.990 

.979 

3.9 

2.9 

1.7 

October  1968 

+63. 1; 

-2.6 

-22.7 

+65. U 

-2.5 

-21.1; 

18.6 

12.9 

10.2 

18.7 

13.8 

9.6 

10. U 

5.8 

3.0 

.996 

.989 

.975 

3.2 

2.0 

1.3 

February  1 969 

+80.0 

-0.81; 

-22.1 

+77.8 

-0.98 

-20.1 

20.8 

17.2 

13.2 

20.5 

16. U 

11.7 

11.1 

7.5 

li.7 

.991; 

.987 

.972 

3.7 

2.8 

2.0 
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6.  VERIFICATION  OF  TECHNIQUE 

The  zero  line  values  were  verified  for  February  1 968  over  the  area 
shown  in  figure  6.  For  each  day  the  position  of  the  zero  line  was 
inferred  from  the  cloud  photographs.   Latitude-longitude  locations  along 
the  zero  line  were  noted  at  half-gridlength  (NMC  grid  scale) .  At  these 
locations  values  of  the  first-guess  long-wave  component  field  were  com- 
pared to  the  total  height  indicated  by  the  operational  analysis.  The 
same  comparisons  were  made  at  locations  taken  from  a  "random"  zero  line. 
This  line  was  drawn  to  represent  the  synoptic  climatology  of  the  area. 
Finally,  at  the  locations  taken  from  both  the  "inferred"  and  "random"  zero 
lines  the  total  height  of  the  first-guess  was  compared  to  the  final 
analysis.  The  results  are  presented  in  table  U. 

Table  U. — Zero  line  verification  for  February  1968 


"Inferred"  zero  line — 
"Random"  zero  line 


RMS  discrepancy 
between  inferred 
zero  line  heights 
and  final  analysis 


RMS  discrepancy 
between  first- 
guess  heights 
at  inferred  zero 
line  locations 
and  final  analy- 
sis 


Average  gra- 
dient of 
final  analy- 
sis at 
inferred 
zero  line 
locations 


Meters 
73 


Meters 
31 
33 


Meters  per 
gridlength 
107 
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The  left-hand  column  of  table  h   indicates  the  usefulness  of  the  cloud 
photographs  in  determining  the  location  of  the  zero  line.   The  "inferred" 
zero  line  heights  represent  a  considerable  improvement  over  the  "random" 
zero  line  heights.  A  more  important  question,  however,  is  the  accuracy 
of  the  "inferred"  heights  derived  along  the  "inferred"  zero  line.   Com- 
paring the  left  and  middle  columns  suggests  that  "inferred"  heights  are 
less  accurate  than  the  coincident  values  of  the  first-guess  field,  but 
the  conclusion  would  be  unwarranted.  The  National  Meteorological  Center's 
(MC)  objective  500-mb  analysis,  used  as  the  verification  standard,  is 
biased  toward  the  first-guess  in  inverse  proportion  to  the  amount  of 
available  data;  and  the  area  chosen  for  verification  is  not  a  dense  data 
region.  It  is  also  apparent  from  the  right-hand  columm  in  table  h   that 
the  zero  line  occurs  in  areas  of  very  tight  gradient,  and  shifting  a 
contour  in  the  total  height  field  by  only  half  a  gridlength  is  sufficient 
to  explain  the  average  discrepancy  between  the  "inferred"  heights  and  the 
final  analysis.  On  this  basis  we  believe  the  indicated  difference  between 
the  "inferred"  and  analyzed  heights  to  be  within  an  acceptable  range  of 
uncertainty. 

As  shown  in  figure  10  and  table  ^,   the  magnitude  of  the  discrepancy 
between  zero  line  and  analyzed  height  values  is  dependent  on  location 
relative  to  the  cloud  system.  The  region  of  largest  discrepancies  occurs 
in  the  main  cloud  shield.  We  found  this  to  be  caused  by  basic  incompati- 
bility between  the  scale  of  our  interpretation  of  the  zero  line  and  the 
resolution  of  the  analysis  system.  The  tight  curvature  of  the  zero  line 
suggested  by  the  cloud  patterns  was  not  resolved  on  the  grid  used  for 
objective  analysis. 

Table  5. — Zero  line  verification  according  to  position  within  cloud 
system  for  "inferred"  heights  differing  from  the  final 
analysis  by  at  least  30  meters 


Ratio  of  estimates  RMS  error  of 
Position  with  discrepancy    estimates  with 

of  30  m  or  more     discrepancy  of 
to  total  sample 30  m  or  more 


Meters 


Main  cloud  shield 

Cold  frontal  band 

Secondary  development 

Boundary  of  open-celled  convection 


15U/2142 

83 

52/290 

60 

117/281 

62 

282/1181 

58 
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Figure  10. --Areas  within  cloud  systems  where  the  zero  line  heights 
were  verified:   Area  I,  main  cloud  shield^  Area  II,  cold 
frontal  band;  Area  III,  secondary  development;  Area  IV,  boundary 
of  open-celled  convection. 
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Figure   11. --Scatter  diagram  of  the  actual  versus   the  residual  of 
the  intensity  of  the  short-wavelength  system. 
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The  area  where  the  second  largest  relative  niunber  of  discrepancies 
occurs  is  in  the  region  of  secondary  systems,  such  as  vorticity  centers 
in  cold  air  or  vorticity  centers  associated  with  enhancement  of  main 
frontal  cloud  bands  (see  "Whitney  and  Herman,  1969).   In  these  situations 
the  accuracy  of  both  the  first-guess  and  final  analysis  fields  is  most 
suspect  in  sparse-data  regions.   Consequently,  we  feel  that  the  zero  line 
data  are  making  their  greatest  contribution  in  these  regions,  even  though 
the  apparent  error  is  large. 

The  regression  equation  for  determining  the  negative  extrema  was 
verified  with  an  independent  set  of  data  consisting  of  131  cases  from 
February,  September,  October,  and  November  1969.  Results  with  the 
independent  sample  are  compared  to  the  dependent  sample  in  table  6. 


Table  6. — Regression  equation  and  sample  statistics  for  determining 
negative  extrema  of  the  short-wavelength  field  (meters) 


Sample 

Maximum 

Minimum 

Mean 

Standard 
deviation 

Regression 
error  of 
estimate 

Regression 

correlation 

coefficient 

Dependent 

Independent— 

-20 

-5 

-322 
-2l|8 

-139 
-139 

61 
58 

38 

.80 
.72 

Table  6  shows  that  the  regression  equation  is  stable.   The  accuracy  of  the 
heights  derived  from  the  regression  equation  is,  as  with  the  zero  line 
heights,  difficult  to  verify.  The  regression  equation  does  produce  some 
systematical  error.  A  scatter  diagram  of  the  predicted  minus  the  actual 
values  versus  the  actual  values  for  the  independent  sainple  is  shown  in 
figure  11.  This  diagram  shows  a  consistent  bias  in  the  predicted  values: 
the  cases  that  are  more  intense  than  the  mean  are  under-predicted,  whereas 
those  weaker  than  the  mean  are  over-predicted.  This  appears  to  be 
characteristic  of  a  least-squares  fit  to  a  sample  of  data,  and  efforts  to 
correct  for  this  bias  have  not  shown  appreciable  success. 

Appendix  A  provides  the  details  adapting  these  techniques  to  a  real- 
time operation  and  illustrates  the  outcome  with  three  cases. 

7 .   CONCLUSIONS 


It  has  been  demonstrated  that  quantitative  use  of  satellite-observed 
cloud  patterns  can  be  made  in  modifying  the  500-mb  height  field.  The 
cloud  pictures  permit  the  identification  and  location  of  the  short- 
wavelength  systems  in  the  500-mb  height  field.  Heights  along  the  zero 
line(s)  in  the  relative  vorticity  field  are  easily  derived  and  demonstrate 
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reasonable  accuracy.  The  regression  equations  permit  quantitative  esti- 
mates of  the  magnitudes  of  the  negative  extrema  in  the  short-wavelength 
field,  and  conversion  of  these  values  to  absolute  heights  gives  reasonable 
estimates  of  the  500-mb  heights  at  the  location  of  the  spiral  cloud 
centers.  It  is  concluded  that  the  technique  can  significantly  improve 
the  accuracy  of  the  500-mb  height  fields  in  sparse-data  areas. 
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APPENDIX  A 

Application 

At  present,  these  zero-line  and  regression  techniques  are  being 
applied  only  to  the  0000  GMT  analysis  in  the  area  covering  the  North 
Pacific  Ocean  from  about  the  International  Date  Line  to  the  west  coast 
of  North  America.   The  12-hour  numerical  ^OO-mh   height  forecast  from  the 
previous  1 200  GMT  analysis  is  computer  processed  to  produce  the  short- 
and  long-wavelength  fields.   These  products  are  available  to  the  analyst 
at  about  2200  GMT.   The  primary  satellite  data  used  are  Automatic  Picture 
Transmission  System  (APT)  pictures  of  cloud  patterns  over  the  eastern  and 
central  Pacific  Ocean.   These  data  received  from  the  satellite  at  San 
Francisco  and  Hawaii,  are  relayed  directly  to  the  National  Environmental 
Satellite  Center  (NESC),  Suitland,  Md.   Although  these  data  provide 
adequate  coverage  to  perform  the  analyses,  they  are  from  3  to  6  hours 
earlier  than  the  0000  GMT  analysis  time.   This  is  partially  compensated 
for  by  use  of  2300  GMT  data  from  the  geosynchronous  Applications 
Technology  Satellite,  ATS  1 .   A  preliminary  analysis  is  performed  using 
the  APT  data  and  then  adjusted  by  using  the  ATS  1  observation. 

The  500-mb  heights  at  the  spiral  cloud  centers  and  along  the  zero 
line(s)  are  derived  and  plotted  on  a  1:20,000,000  polar  stereographic 
map.   The  locations  of  troughs,  ridges,  and  the  axes  of  maximum  anti- 
cyclonic  shear  inferred  from  the  cloud  patterns  also  are  entered  on  this 
map.  Using  the  derived  500-mb  heights,  the  synoptic  features  indicated 
above,  and  the  general  configuration  of  the  cloud  patterns,  the  analyst 
constructs  a  contour  height  analysis  for  the  area. 

The  completed  analysis  is  given  to  'the  Automated  Analysis  Branch  of 
NMC  at  approximately  0100  GMT;  these  data  are  then  blended  with  the 
0000  GMT  radiosonde  observations  through  NMC's  ^OO-mb  hand  analysis.   The 
final  hand  analysis  is  compared  to  a  preliminary  machine  analysis  (the 
RADAT  analysis),  and  NMC  analysts  generate  "bogus"  wind  and  height  reports 
to  correct  the  discrepancies.   These  "bogus"  heights  are  entered  into  the 
operational  machine  analysis  which  is  used  as  the  initial  data  for  the 
Primitive  Equation  Model  forecast. 

Satellite  Input  to  Numerical  Analysis  and  Prediction  (SINAP)  has  been 
used  daily  since  July  1,  1969.   This  section  illustrates  several  examples 
of  the  operational  product. 

0000  GMT  October  16,  1 969 

Figure  12,  top,  shows  the  ESSA  9  digitized  mosaics  for  approximately 
0000  gmt,  October  I6,  I969.   Superimposed  on  the  mosaic  are  the  first- 
guess  field,  500-mb  contours,  and  the  zero  line  of  the  short-wavelength 
field  as  derived  from  the  first-guess  field.   The  discrepancies  between 
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Digitized  cloud  mosaic,  500-mb  first-guess  contours  (decameters)  and 
short-wavelength  field  zero  line. 


SINAP  analysis 

Figure  12. --Cloud  pattern,  numerical  fields,  and  SINAP  analysis  for 
0000  GMT,  October  16,  I969. 
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the  zero  line  and  the  cloud  patterns  shown  in  the  figure  are  rather 
typical.   The  general  pattern  of  the  zero  line  reflects  the  systems  which 
appear  in  the  cloud  patterns,  but  there  are  errors  in  location  and  ampli- 
tude.  For  example,  the  ridge  line  near  170°  W.  appears  too  far  east 
compared  with  the  cloud  patterns.   Near  i|0°  N.  and  l65°  E.  the  zero  line 
incorrectly  extends  through  an  area  of  inversion-dominated  stratocumulus 
clouds.   Near  U0°  N. ,  1^0°  W. ,  the  zero  line  is  incorrectly  located  in  the 
middle  of  a  field  of  open-celled  convective  cloudiness,  the  SINAP  analysis 
for  0000  GMT,  October  16,  I969  is  shown  in  figure  12,  bottom.   This 
analysis  contains  the  heights  along  the  zero  line  (in  decameters  with  the 
first  digit  deleted)  and  at  three  vorticity  centers  computed  from  the 
regression  equation  (decameters).   The  height  at  the  vorticity  center  located 
near  k.k°   N.,  179°  W.  was  5jUlO  meters;  this  value  differed  from  the  first- 
guess  field  by  -l50  meters.   Heights  at  the  two  vorticity  centers  located 
near  ship  PAPA  (50°  N.,  Ili5°  W. )  were  5,220  meters  and  5:, 270  meters.   These 
heights  deviated  from  the  first-guess  field  by  about  -60  meters. 

In  figure  13?  top,  the  500-mb  contours  from  NMC's  hand  analysis  are 
superimposed  upon  the  ESSA  9  digitized  mosaic.   Also  shown  is  the  zero 
relative  vorticity  line  derived  from  the  machine  analysis  for  0000  GMT, 
October  16, I969.   A  reasonable  correspondence  between  the  cloud  patterns 
and  the  zero  relative  vorticity  line  is  apparent,  although  improvements 
could  still  be  made.   (Note  that  the  zero  line  still  encloses  a  large 
area  of  stratocumulus  near  IlO°  N.,  175°  E.) 

In  figure  13,  bottom,  the  SINAP  heights  are  shown  together  with  the 
available  500-mb  wind  and  height  reports.  Note  that  the  zero  line  heights 
near  radiosonde  stations  appear  quite  compatible  with  the  reported  heights. 
In  particular,  the  location  of  the  zero  line  and  its  derived  heights  is 
completely  compatible  with  the  apparent  location  of  the  jet  stream  as 
indicated  by  the  wind  speed  at  the  station  on  the  east  coast  of  Kamchatka 
Peninsula.  No  data  are  available  to  verify  the  system  located  near  Ulj.°  N. , 
179°  ¥.;  the  height  at  the  vorticity  center  located  near  50°  N. ,  l50°  W. , 
however,  was  compatible  with  the  observed  height  and  wind  at  ship  PAPA. 

0000  GMT,  February  9,    1970 

The  previous  case  represents  a  classical  example  where  the  distinct 
change  between  cloud  types  leaves  little  doubt  as  to  the  proper  location 
of  the  vorticity  features.   This  is  usually  the  exception  rather  than 
the  rule;  consequently,  the  analyst  must  make  subjective  decisions  as  to 
the  locations  of  the  zero  line  and  the  vorticity  centers .   This  example 
and  the  next  were  chosen  at  random  from  the  daily  operational  analyses. 
The  complex  cloud  patterns  of  these  cases  represent  the  usual  conditions 
which  confront  the  analyst.   The  SINAP  analyses  shown  have  been  slightly 
modified  from  the  operational  product  primarily  for  esthetic  purposes. 
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Digitized  cloud  mosaic,  NMC's  operational  500-mb  analysis  and 
500-mb  relative  vorticity  zero  line 


SINAP  heights  and  radiosonde  data 

Figure  13. --NMC's  operational  analysis,  SINAP  heights,  and  radiosonde 
data  for  0000  GMT,  October  16,  1969. 
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but  also  to  insiire  registration  with  the  digitized  ATS-1   photograph^; 
they  differ  only  in  minor  details  from  the  operational  products  which 
were  produced  independently  of  the  radiosonde  data. 

Figure  ih,    top,  shows  the  first-guess  500-mb  height  contours  and  the 
zero  line  in  the  short-wavelength  field  for  0000  GMT,  February  9,  1970, 
superiinposed  on  the  2130  GMT,  February  8,  1970,  ATS-1  digital  product. 
The  first-guess  field  in  this  case  shows  good  agreement  with  the  large- 
scale  pattern  of  the  cloudiness,  but  contains  no  information  on  the  short- 
wavelength  systems  indicated  by  the  cloud  patterns.   Conspicuous  by  its 
absence  from  the  contour  field  is  the  vorticity  maximum  indicated  by  the 
intense  vortical  cloud  pattern  located  near  14.2°  M. ,  l60°  W.  Another 
deficiency  in  the  first-guess  field  is  the  trough  which  connects  the  deep 
low  centered  near  50°  N.,  175°  W.  to  the  subtropical  system  located  near 
35°  N.,  135°  W.  The  presence  of  this  trough  places  the  area  N.  liO°-50° 
and  W.  135°-l55°  under  a  positive  vorticity  regime,  whereas  the  cloud 
patterns  indicate  the  opposite. 

The  SINAP  analysis  is  shown  in  figure  ^h,   bottom.  Three  vorticity 
maxima  were  entered  into  the  SINAP  analysis  which  were  not  present  in 
the  first-guess  field.  These  systems  were  located  near  1|0°  N. ,  175°  W.; 
U2°  N.,  160°  W.;  and  26°  N.,  130°  W.   The  heights  computed  at  these 
centers  deviated  from  the  first-guess  field  by  -60,  -60,  and  +lj.O  meters 
respectively.   In  the  SINAP  analysis  the  pattern  of  the  positive  vorticity 
was  separated  into  two  separate  pools;  one  associated  with  the  deep  low 
in  the  Central  Pacific  and  the  second  with  the  cut-off  low  west  of 
California.   Strong  ridging  and  negative  vorticity  were  indicated  in  the 
region  where  the  first-guess  field  indicated  a  trough.  Another  area  where 
considerable  adjustment  was  indicated  was  near  30°  N. ,  175°  W.;  here  the 
first-guess  field  zero  line  was  adjusted  some  5  degrees  to  the  south. 

The  operational  500-mb  height  analysis  and  the  zero  relative  vorticity 
line  derived  from  it  are  shown  superimposed  upon  the  ATS-1  digitized 
pictures  in  figure  l5^  top.  It  is  clear  that  many  details  of  the  SINAP 
analysis  have  been  incorporated  into  NMC  numerical  product.   A  secondary 
center  has  been  entered  into  the  analysis  in  the  southeast  quadrant  of  the 
cut-off  low  in  the  Eastern  Pacific.   A  sharp  short-wavelength  system  has 
been  analyzed  near  i|5°  N.,  155°  W.  which  corresponds  reasonably  with  the 
cloud  patterns .  Of  primary  importance  is  the  ridge  which  has  been  entered 
into  the  analysis  along  ll4.5°  W. ;  this  has  effected  a  reversal  in  the  sign 
of  the  vorticity  of  the  first-guess  field  in  this  area.  The  zero  line 
between  150°  W.  and  160°  W.  shows  some  noise;  but  with  gradients  of  the 
magnitude  indicated  by  the  analysis  in  this  area,  very  subtle  changes  in 


^This  picture  is  produced  daily  from  the  2130  GMT  ATS-1  cloud  photograph. 
The  cloud  photograph  is  processed  by  machine  and  remapped  in  full  resolution 
on  a  1:20,000,000  polar  stereographic  projection. 
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Digitized  cloud  photograph,  500-mb  first-guess  contours  (decameters) 
and  short-wavelength  field  zero  line 


SINAP  analysis 

Figure  1^. — Cloud  pattern,  numerical  fields,  and  SINAP  analysis 
for  0000  GMT,  February  9,  1970. 


Digitized  cloud  photograph,  NMC's  operational  500-mb  analysis  and 
500-mb  relative  vorticity  zero  line 


\^ 


SINAP  height  and  radiosonde  data 

Figure  15 . — MC ' s  operational  analysis,  SINAP  heights,  and  radiosonde 
data  for  0000  GMT,  February  9,  1970. 
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the  curvature  of  the  contours  can  produce  large  changes  in  the  cor- 
responding vorticity  field.  This  noise  is  considered  of  negligible 
significance. 

The  SINAP  heights  and  the  actual  radiosonde  and  reconnaissance  data 
are  shown  in  figure  ]$,   bottom.  With  the  exception  of  one  area,  the  heights 
along  the  zero  line  and  at  the  vorticity  centers  appear  perfectly  com- 
patible with  the  reported  heights.  The  SINAP  heights  along  the  trailing 
edge  of  the  cloud  band  associated  with  the  cut-off  low  are  consistently  too 
high.  This  is  a  factor  that  has  been  observed  consistently  in  the  presence 
of  intense  convective  cloud  bands  associated  with  this  type  of  system. 
Placing  the  zero  line  along  the  leading  edge  of  such  a  cloud  band  produces 
heights  which  are  much  more  compatible  with  observations. 

0000  GMT,  March  23,    1970 

The  500-mb  first-guess  height  field  and  the  zero  line  in  the  short- 
wavelength  field  are  shown  superimposed  upon  the  ATS-1  digitized  product 
for  2130  GMT,  March  22,  1970  in  figure  16,  top.  As  in  the  previous  cases, 
the  contour  pattern  in  the  first-guess  field  and  the  location  of  the  zero 
line  correspond  reasonably  with  the  cloud  patterns.  In  this  case  the 
major  deficiency  in  the  first-guess  field  is  the  lack  of  amplitude  in  the 
short-wavelength  systems. 

The  satellite  photographs  at  this  time  showed  an  intense  "comma" 
shaped  cloud  pattern  at  ilu°  N.,  179°  E.  The  SINAP  analysis  in  figure  l6, 
bottom,  correctly  included  this  cloud  pattern  entirely  within  the  area 
of  positive  vorticity.   This  bright  comma-shaped  cloud  pattern  is  associated 
with  an  area  of  strong  positive  vorticity  advection  within  the  cold  air  in 
advance  of  the  500-mb  trough  line.   The  northward  bulge  of  the  frontal  band 
in  the  vicinity  of  h$°   N. ,  l60°  W. ,  indicates  the  vorticity  maximum  is 
interacting  with  the  primary  baroclinic  zone.   In  such  transition  situations 
the  location  of  the  zero  line  in  the  area  of  interaction  is  not  always 
distinct,  but  the  correct  analysis  is  to  place  the  zero  line  along  the 
northwest  boundary  of  the  dense  cirriform  cloudiness  associated  with  the 
frontal  band.  The  sag  in  the  contours  located  near  33°  N. ,  175°  W.  in  the 
first-guess  field  produced  a  marked  disagreement  between  the  cloud  patterns 
and  the  first-guess  field  zero  line.   The  clouds  show  that  this  feature  was 
erroneous;  it  was  eliminated  in  the  SINAP  analysis. 

NMC's  operational  500-mb  height  analysis  for  0000  GMT,  March  23,  1970 
and  the  corresponding  zero  relative  vorticity  line  are  shown  superimposed 
upon  the  2130  GMT,  March  22,  1970,  ATS-1  digitized  picture  (figure  17,  top). 
Comparing  the  contour  patterns  of  figures  1 6  and  1 7  reveals  that  the 
operational  analysis  contains  much  more  detail  than  was  evident  in  the 
first-guess  field.  The  enhancement  of  the  gradients  in  the  operational 
analysis  is  particularly  apparent  compared  with  the  first-guess  field. 
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Digitized  cloud  photography,  500-inb  first-guess  contours  (decameters) 
and  short  wavelength  field  zero  line 


SIMP  analysis 

Figure  16. — Cloud  patterns,  niJimerical  fields,  and  SINAP  analysis  for 
0000  GMT,  March  23,  1970. 
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Digitized  cloud  photograph,  NMC's  operational  ^00-mb  analysis  and 
500-mb  relative  vorticity  zero  line 


SINAP  heights  and  radiosonde  data 

Figure  1?. — NMC's  operational  analysis,  SINAP  heights,  and  radiosonde 
data  for  0000  GMT,  March  23,  1970. 
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The  zero  relative  vorticity  line  as  derived  from  the  operational  height 
analysis  is  very  noisy,  however;  and  the  analysis  still  contains  the 
erroneous  area  of  positive  vorticity  centered  near  35°  N. ,  172°  W.   This 
is  attributed  to  an  excess  of  cyclonic  curvature  of  the  contours  in  this 
area. 

In  figure  17,   bottom,  are  the  heights  derived  from  the  SINAP  analysis, 
together  with  the  available  radiosonde  data  for  0000  GMT,  March  23,  1970. 
There  are  only  three  points  where  the  SINAP  data  were  sufficiently  close 
to  radiosonde  reports  to  allow  direct  comparison  of  the  heights.   The 
height  along  the  zero  line  2  degrees  east  of  ship  PAPA  was  5jU60  meters; 
the  reported  height  at  PAPA  was  ^ ,h90   meters,  and  the  wind  was  260°  at 
85  knots.  A  height  of  5^830  meters  was  derived  along  the  zero  line  some 
2°   east  of  Hawaii;  this  was  10  meters  higher  than  the  reported  height  at 
Hilo.   A  height  of  5,700  meters  was  computed  at  the  open  vortex  located 
at  27°  N. ,  lU2°  W. ;  this  appears  consistent  with  the  height  reported  by 
ship  NOVEMBER  (30°  N. ,  1^0°  W) . 
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APPENDIX  B 
Properties  of  Scale  Separation 
The  smoothing  filter  is  defined  by 


0  =  4>0  +  8^J^''   vVda  (B1) 

where  <f>  is  the  smoothed  field,  00^^  '^^^  original  field,  5  is  a  constant, 
and  a  is  the  degree  of  smoothing  parameter.  For  purposes  of  later  gener- 
alization it  is  convenient  to  express  equation  (B1 )  in  operator  notation: 


<p[a]  =  0[o]-0  [0,a] 


(B2) 


where  0 [a]  is  to  be  read  as  "the  field  <p  smoothed  to  degree  of  smoothing 

a^"  and0[o,ajis  to  be  read  as  "the  residual  field  resulting  from 
smoothing  the  <f>   field  from  the  state  a  =  0        to  a  =  a  . " 

In  order  to  investigate  the  properties  of  the  smoother,  let  us 
assume  that  the  field  (p  is  represented  by  a  one -dimensional  wave  with 
wave  number  k.  The  amplitude  of  the  wave  will  be  a  function  of  the  degree 
of  smoothing.  In  differential  form  equation  (B1 )  is  given  by 


ba  ^  (B3) 

and  the  one-dimensional  wave  is  expressed  by 


^[o]  =  A,,  [ole*"^-  (bU) 
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Substituting   (BH)   into    (B3)   and  solving  gives 


0[a]=0[o]e-^'5'« 
Ak[a]  =  Ak  [O]  e 


w'8'a 


The  amplitude  reduction  factor  caused  by  the  smoothing  is  given  by 
Ak[0] 


(B5) 


(B6) 


(B7) 


It  is  readily  seen  that  the  amplitude  reduction  is  a  function  of  the 
degree  of  smoothing  and  the  wave  number.  The  shorter  the  wave  (large  k) 
or  the  greater  the  smoothing  (larger  a)   the  more  the  amplitude  will  be 
reduced. 

In  practice  equation  (B1 )  is  integrated  by  repeated  applications  of 
the  del  square  operator  with  varying  increments  in  a   (ai,a2,-".^n)  •   It,  can 
be  seen  from  equation  (B5)  that  the  resultant  smoothing  is  cumulative  in 
a     since 


0[a,][a,]  =  0[o][aJe-'<'«'<''. 

In  explicit  finite  difference  form: 


(B8) 


0[a,]  =0[O]+V'0[o]a, 
0[a,+  a2]  =  0[a,]-t-v2  0[a,]a; 


=  i> 


-n-i     -1 


0 


n-1 


Li-i       . 


(B9) 
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where  'iP'^  is  the  five-point  finite  difference  del  operator  on  a  lattice 
with  gridlength  S  .  In  general  the  a^  are  not  equal,  but  for  reasons 
of  numerical  stability  increase  with  successive  steps. 

Smoothing  is  complete  when  a=S^i   reaches  a  predetermined  value. 
0  [a]  is  then  the  long-wave  component  field,  and  the  short-wavelength 
component  field  is  formed  by  subtracting  the  long-wave  from  the  original 
field. 

The  degree  of  smoothing  is  assigned  by  deciding  on  the  reduction 
factor  to  be  applied  to  a  particular  wavelength.   In  order  to  elucidate 
this  point,  it  is  convenient  to  construct  the  general  smoothing  curve 
(filter)  from  equation  (B7)  with  the  wave  number  replaced  by  the  wave- 
length in  grid  increments. 


K  =27T/n8  (BIO) 

where  n  is  the  number  of  gridlengths  corresponding  to  K.   Then 


R  =  e-^^'«/"'   "  (BID 

Using  equation  (B11)  it  is  possible  to  plot  a  series  of  curves  cor- 
responding to  different  n  for  the  reduction  factor  as  a  function  of  a  . 
However,  a  single  general  curve  can  be  constructed  by  introducing  a  new 
variable  x  which  is  a  function  of  both  a  and  n: 


rx^  =  47r^a/n2  (B12) 

where  r  is  an  arbitrary  positive  constant  used  to  scale  x  (as  will  become 
apparent) . 

Thus  from  (fill) 


R  =  e-^^^  (B13) 
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Figure  l8. — Normalized  filter  curves  for  first-  and  third-order 
smoothing. 
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For  plotting  purposes  it  is  convenient  to  define  r  by  requiring  that 

R  =  .5  when  x  =  1.  Using  this  definition  in  equation  (B^3),    r  =  .6931^7. 

The  filter  curve  is  depicted  by  the  solid  line  in  figure  18. 

To  demonstrate  the  use  of  the  ciorve,  suppose  that  we  decide  to  apply 
the  smoother  until  a  10-gridlength  wave  is  reduced  to  10  percent  of  its 
original  amplitude.   From  figure  l8  we  find  that  for  R  =  .1,  x  =  1.82. 
Using  equation  CB12)  and  the  definition  of  r,   we  find  that  the  proper 
degree  of  smoothing  .is  o.  =   5.8. 

Because  n  is  linearly  dependent  on  x,  the  slope  of  the  curve  in  figure 
l8  depicts  the  sharpness  of  the  filter.  A  perfect  filter  would  be  a 
vertical  line  at  x  =  1 .  All  shorter  wavelengths  (smaller  n,  larger  x) 
would  be  entirely  removed  (R  =  O),  whereas  all  longer  wavelengths  would  be 
entirely  retained  (R  =  1 ) .   The  filter  described  by  the  smoothing  con- 
sidered above  (hereafter  called  first-order  smoothing)  is  not  particularly 
sharp.   The  fact  that  the  filter  is  not  sharp  means  that  the  scale 
separation  is  imperfect.  The  long -wavelength  field  contains  some  com- 
ponents of  the  short-wavelength  field  and  vice  versa.   This  undesirable 
feature  can  be  corrected,  at  the  expense  of  computation  time,  by  applying 
higher  order  filters . 

The  principle  of  higher  order  smoothing  is  conceptually  simple.   For 
example,  suppose  that  the  short-wavelength  field  after  first-order 
smoothing  contains  95  percent  of  the  amplitude  of  a  particular  short  wave 
and  erroneously  contains  15  percent  of  the  amplitude  of  a  particular  long 
wave  (present  because  of  the  imperfection  of  the  filter  associated  with 
first-order  smoothing).   If  after  completion  of  first-order  smoothing  the 
short-wavelength  field  is  further  smoothed  into  large  and  small  scale 
components,  the  amplitude  of  the  long  wave  will  be  reduced  in  the  second 
short-wavelength  component  field  to  .15  X  .15  or  2.25  percent  of  its 
original  amplitude.  The  amplitude  of  a  short  wave  will  be  reduced  in  the 
second  short-wavelength  component  field  to  .  25  X  .95  or  90.25  percent  of 
its  original  amplitude.   This  second-order  smoother  is  certainly  sharper 
than  the  first-order  smoother  because  the  longer  wavelengths  are  filtered 
more  rapidly  than  the  shorter.   In  fact,  as  will  now  be  shown  analytically, 
the  repeated  smoothing  of  the  residual  field  can  be  extended  to  give  a 
filter  of  any  desired  sharpness. 

Equation  (B2)  for  first-order  smoothing  can  be  rearranged  to  give 


0[O,a>0[o]-0[a]  (Bll;) 
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Second-order  smoothing  as  described  above  involves  smoothing  the  residual 
field  a  second  time: 

0[O.a][O,a]=  i<A[o]-0[a]|[O,a] 

=  0[O][o]-0[o][a]-<^[a][o]  +  0[a][a]  (Bl5) 

From  equation  (B8) 


0[O,a]rO,a]=0[o]-20[a]  +  0[2a]  (Bl6) 

Adopting  the  convention 


0  [a]  [a]  =  0  [2  a]  =(/>[«]'  (B17) 

equation   (Bl6)    can  be  written 


(B18) 


0[O,af  =0(l-[a])' 

or  in  general 

0[O,a]''=0(l-[a])'' 


The  final  long  wave  component  field  after  nth  degree  smoothing  is   then 
found  by  subtracting  the  residual  from  the  original  field: 


•i^Lw^  ^  [O]-0[O.  «]"=  <?^     l-(l-[a])''  (B19) 
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For  reasons  of  economy  we  have  adopted  third-order  smoothing.   The  cor- 
responding long-wavelength  field  is  readily  shown  to  be 


0^^  =  30[a]-30[2a]+9!)[3a]  (B20) 

As  shown  by  equation  (B20),  the  long -wave length  field  is  found  in  practice 
not  by  repeatedly  smoothing  the  residual  field  but  by  smoothing  the 
original  field  to  3  at  .   This  is  computationally  more  efficient  because, 
as  mentioned  earlier,  as  smoothing  progresses  larger  smoothing  increments 
can  be  taken  (viz.,  equation  (B9)). 

Using  equations  (B5),  Cb6),  (B7),  (BIO),  and  (B12)  with  (B20),  the 
reduction  factor  for  third-order  smoothing  is 


V  R  =  3e-'''''-3e-2'-xVe-3''x2  (g2i ) 

Again  adopting  the  convention  that  R  =  .5  when  x  =  1,  r  =  1  .578I|.261|.   The 
third-order  filter  curve  is  plotted  as  a  dashed  line  on  figure  l8.   As 
anticipated,  the  third-order  filter  is  considerably  sharper  than  the  first- 
order  filter. 

To  use  the  third-order  smoother,  suppose  as  before  that  we  wish  to 
reduce  a  10-gridlength  wave  to  10  percent  of  its  original  amplitude.   From 
the  figure  we  find  that  from  R  =  .1,  x  =  1  .1|6.   Using  equation  (B12)  and 
the  definition  of  r,  we  find  that  the  proper  degree  of  smoothing  is 
a   =   8.52. 

Equation  (B1 )  describes  the  procedure  for  explicit  smoothing.   In 
practice  ijnplicit  smoothing  is  found  to  be  more  efficient.   One  smoothing 
cycle  is  represented  by 


The  field  is  solved  by  standard  overrelaxation  using  a  relaxation  factor 


(Holl,  1967)  (o   equal  to 


CJ  =  ]  4 


U(l-q2) 


V7 


where 


(B23) 


q  =  25a/(25a+l) 


(B2U) 


The  8a    for  each  smoothing  cycle  is  assigned  such  that  15  cycles  will  give 
fields  smoothed  to  a  ,  2  a  ,  and  3  a  •   A  normalized  table  for  5a*   is 
given  in  table  7.   The  actual  80.     used  in  each  cycle  is  found  by 
multiplying  the  normalized  value  by  the  desired  degree  of  smoothing. 

Table  7. — Normalized  smoothing  increments  (from.  Holl,  I967) 


Cycle  No. 

5  a* 

la* 

Cycle  No. 

5a* 

25a* 

1 

.006 

.006 

9 

.15U 

.Uii8 

2 

.009 

.015 

10 

.230 

.678 

3 

.Olli 

.029 

11 

.322 

1.000 

h 

.020 

.oa9 

12 

.500 

1.500 

5 

.030 

.079 

13 

.^00 

2.000 

6 

.Oli5 

.12U 

1U 

.500 

2:500 

7 

.068 

.192 

15 

.500 

3.000 

8 

.102 

.29U 

— 

— 

—  — 
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APPENDIX  C 
Program 

SUBROUTINE  SCALE 

DIMENSION  Fl(302l),F2(302l),F3(302l),MFAC(302l),ZOUT(302l) 

COMMON  Fl,F2jF3,MFAC,Z0UT 

REAL  MFAC 
C  IMPLICITE  SCALE  SEPARATION  PROGRAM  DESIGNED  BY  M.M.HOLL 

C  PROGRAMMED  BY  C.HAYDEN  AND  R.NAGLE  NOAA  NESS 

C  FIELD  TO  BE  SMOOTHED  ENTERS  INTO  Fl( ORIGINAL)  AJ^D  IS  OPERATED  ON 

C  BY  SUBROUTINE  SCLN(N,II)   FROM  F2.    ORIGINAL  ETELD  MUST  BE  A  53  XH7 

C  ARRAY  AND  SCALED  IN  METERS.    NINETY-FIVE (95)   PERCENT  OF  THE 

C  AMPLITUDE  OF  ¥A¥E  NUMBER  10   (OR  HIGHER)  IS  ELIMINATED  FROM  THE 

C  ORIGINAL  FIELD  BY  THIS   PROGRAM,   THE  RESULTING  FIELD  IS  CALLED 

C  THE  LONG  WAVE  PATTERN   (SR).   THE  SR  ZTELD   IS  CONTAINED  IN  F3. 

THE  PORTION  OF  THE  ORIGINAL  FIELD  WHICH  HAS   BEEN  REDUCED  BY  THE 
C  SMOOTHING  PROCESS  MAY  BE  RETRIEVED  BY  SUBTRACTING  THE  SR  FIELD 

C  FROM  THE  ORIGINAL,   THIS    (SD)   IS  FOUND  IN  F2. 

C  THE  SMOOTHING  IS  ACCOMPLISHED  BY  OVER-RELAXATION.   THE  ORIGINAL 

FIELD  IS  SMOOTHED  TO  ALPHA=1,   THEN  TO  ALPHA=2  AND  ALPHA=3. 
C  THE  SR  FIELD  IS  PRODUCED  AS  FOLLOWS 

C  SR=3-"-Z(ALPH/l=l)-3-"-Z(ALPHA=2)+Z(ALPHA=3)   WHERE  Z  (ALPHA =l)   REPRE 

C  SENTS   THE  ORIGINAL  FIELD  SMOOTHED  TO  ALPHA=1  AND  SO  ON 

C  NO  INPUT-OUPUT  ROUTINES  ARE  PROVIDED 

CALL  MFACTR 

N=5 
C  SMOOTH  TO  ALPHA1=5  AND  STORE  RESULTING  FIELD  IN  ZOUT 

DO  10  1=1,11 

CALL  SCLN(N,I) 
10  CONTINUE 

DO  20  1=1,3021 

Z0UT(I)=F2(I) 
20  CONTINUE 
C  SMOOTH   TO  ALPHA2=10,    COMBINE  WITH  3-"-Z(ALPHA=1)AND  STORE  IN  ZOUT 

CALL  SCLN(N,12) 

CALL  SCLN(N,13) 

DO  30  1=1,3021 

Z0UT(I)=3.-'KZ0UT(I)~F2(I) 
30  CONTINUE 
C  SMOOTH  TO  ALPHA3=15 

CALL  SCLN(N,m) 

CALL  SCLN(N,15) 
C  FORM  SR  AND  SD  FIELDS 

DO  I4O  1=1,3021 

F3(I)=Z0UT(I)+F2(I) 

F2(I)=F1(I)-F3(I) 
10  CONTINUE 

RETURN 

END 
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SUBROUTINE  MFACTR 

DIMENSION  F1(3021),F2(3021),F3(3021),MFAG(3021),FLAT(3021), 
1Z0UT(3021) 
COMMON  F1,F2,F3,MFAC,Z0UT,FLAT 
REAL  MFAC 
PIB  2=3. 1^1^926/2. 

RC0N=31. 20^35919 
CKC0N=5. -3. 1^415926/180. 

DO  10  J=l,57  ■         _ 

L=(J-l)-53 

DO  10  1=1,53 

K=L+I 

FJ=J-29 

FI=I-27 

R=SQRT  ( FI-;;-51  +FJ  -x-F  J ) 

FLAT  ( K )  =PIB2-2 . -;;-ATAN(r/RCON) 

if(flat(k).lt.ckcon)5^6 
5  sl=sin(gkg0n) 

GO  TO  7 

6  SL=SIN(FLAT(K)) 

7  S=l. 86602/(1. +SL) 
P1AT(K)=SL 

10  MFAG(K)=S-"-S 
25  GONTINUE 

RETURN 

END 
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SUBROUTINE  SCLN(N,II) 

DIMENSION  F1(3021),F2(3021),F3(3021),MFAC(3021),B(3021),DA(15) 

COMMON  F1,F2^F3:,MFAC 

EQUIVALENCE  If3,B) 

REAL  MFAC 

DATA  (DA=.006,  .009,  .Olii,  .02.  ,03,  .0^5,  .068,  .102,  .l5U, .  230,  .322,  .5^ 

1.5,. 5,. 5) 

DAL=DA(ll)-;rN 

DO  10  J  =  2,56 

L  =  (J-l)-x-.53 

DO  10  I  =  2,52 

M  =  L  +  I 

B(M)  =  F2(M)+DAL-x-MFAC(M)-;KF2(M+1)+F2(M-1)+F2(M+53)+F2(M-53) 
l-B-x-F2(M))/2 
10  CONTINUE 

QB=2.-"-DAL/(  2.--DAL+1. ) 

DAL  =  DAL/2. 

QBS  =  QB-;;-QB 

OM  =  QBS/(1.+SQRTF(1.-QBS)) 

OMEGA  =  i.+0M-;;-OM 

KK  =  0. 
12  KK=KK+1 

if(kk.gt.io)go  to  30 

ICON  =  0 

DO  15  J  =  2,56 

L  =  (j-1)  -;;-  53 

DO  15  I  =  2,52 
M  =  L  +  I 

R  =  (DAL-x-MFAC(M)-;KF2(M+1)+F2(M-1)+F2(M+53)+M-53))  +  B(M))/(1.+I;. 
l-»-DAL-x-MFAC(M)  )-F2(M)" 

RR  ~  R'^'R 

IF(RR.GT.25.)IC0N  =  ICON+1. 

F2(M)  -  F2(M)  +OMEGA-X-R 
15  CONTINUE 

IF  (ICON.GT.O)GO  TO  12 
20  RETURN 
30  PRINT  kO 
kO   FORMAT (-X-  SCALE  NOT  CONVERGING-;:-) 

STOP 

END 
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